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Abstract 

We investigate the existence of stable charged metalUc bubbles using the shell 
correction method. We find that for a given mesoscopic system of n atoms of a given 
metal and q <^ n (positive) elementary charges, a metallic bubble turns out to have 
a lower total energy than a compact spherical cluster, whenever the charge number 
q is larger than a critical charge number qc. For a magic number (n — q) of free 
electrons, the spherical metallic bubble may become stable against fission. 
PACS: 36.40.+d, 36.20.Kd, 31.10.+Z, 21.90.+f, 21.60.Cs 

1 Introduction 

Neutral and charged metal clusters consisting of a few 100 to a few 1000 atoms possi- 
bly containing a limited member (< 10) of positive or negative surplus charges received 
considerable scientific interest since 1984. At that time, when studying the formation of 
tiny alkali clusters out of metal vapour, an enhanced production of clusters with certain 
atomic numbers (n = 8, 20, 40, 58, 92, etc.) was observed [1, 2] and was correlated with 
the appearance of shell closures for the motion of the free electrons i.e. the conduction 
(valence) electrons in a spherically symmetric average potential for the itinerant electrons. 
A free electron feels an average potential which is produced by the background of positive 
ions, on the one hand, and by the other free electrons, on the other. 

The most prominent peaks in the mass yield were shown to be due to shell closures 
in a spherically symmetric potential [3-8], whereas tinyer details of the abundance curve 
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were successfully related to secondary shell effects in axially symmetric deformed [9-11] 
and non-axially symmetric deformed [12] potentials. 

The calculations are usually performed using the Strutinsky shell correction method 
[13] or the more involved self-consistent field approach [14, 15]. 

A large amount of beautiful experimental work [16, 17] has been performed since 
the discovery in 1984. The comparison between the experimental and theoretical work 
(see for instance Ref. 10) is in general satisfactory. The greatest part of the work has 
hitherto been devoted to uncharged clusters. Experiments on charged metallic clusters 
were performed by C. Brechignac et al. Li^"^ (n = nr of atoms, g = nr of elementary charge 
units) in Ref. 18 and for Sb^^ in Ref. 19. It was found that in most cases singly charged 
clusters are stable, i.e. their dissociation is endothermic. Clusters with a positive charge 
g > 1 are observed if the number n of constituant atoms is larger than a critical number nl 
which depends on the system considered. Decay by fission and by evaporation of neutral 
or charged fragments compete with each other, the decay by fission being delayed with 
respect to the decay by evaporation. 

Theoretical studies of the decay were pubhshed by F. Garcias et al. [20] using 
a semi-empirical model for the fission of multiply charged metal clusters and by D. Gross 



In all the theoretical studies of mesoscopic metallic clusters it was assumed that 
the groundstate of the cluster corresponds either to a compact spherical or to a compact 
deformed shape. 

As we show in this paper, for large enough charge q, the state of lowest energy of 
a charged metallic cluster may correspond to a spherical bubble. Within the liquid drop 
approximation, the spherical bubble solutions turn out to be unstable versus fission. This 
has been shown in nuclear physics quite some time ago [22]. In a recent work [23], we 
found that spherical nuclear bubbles may be stabilized by shell effects. In this paper, we 
show that the same is true for charged metal clusters if the number of valence electrons 
corresponds to a closed shell. 

In Section 2, we define the theoretical model and in Section 3 we present the results 
we obtained. In Section 4 contains a short summary and a discussion of open questions. 



For mesoscopic metaUic clusters containing from 100 to more than 1000 atoms, the distri- 
bution of the positively charged ions can be approximately described by a homogeneously 
smeared-out density which in the simplest case is given by a step function. For a bubble 
with inner radius R2 and outer radius Ri, the density of positive ions is thus given by 



[21]. 



2 Theoretical model 



Pion(r) 



ion 



9o{r-R2)9o{Ri-r) , 



[2.1] 
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The constant bulk density Pion> is usually given in terms of the radius of a sphere 
which contains 1 atom on the average 

rs= [3/47rp,jV3. (2.2) 

For a given number n of atoms, the volume of bubble layer has to be equal to the one of 
a compact spherical cluster of radius Rq 

^{Rl - Rl)Kon= Y^o Pion= n . (2.3) 

Consequently, the shape of the spherical bubble is determined by only one free parameter. 
We choose it to be the ratio 

/ := Rl/Rl (2.4) 

between the volume of the inner hole to the volume of the entire bubble. 

The conduction electrons move independently in an average potential V{r) which 
represents the mean interaction of a given electron with all the other electrons and with 
the positive back ground charge Pion('")- In the Hartree-Fock approximation the single 
particle states ^v{r) of the electron^ and the corresponding single particle energies are 
obtained as the selfconsistent solution of the coupled equations 
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'^u{r) = Su^Puir) , (2.5) 



V{r) ■ (^.(f) = [Mon(r) - K(r)] ■ , (2-6) 

VUr) = - [ dr'^^ , (2.7) 



Ve{r) ■ Mr) = E / d''^' fiir')—^ Mr')Mr) - Mr')Mr) ■ (2.8) 

K^U I ' 'I 

In (2.8), Co is the elementary charge and is the occupation number of the single particle 
state For temperature T = 0, is given by 

= 1 for e^^KSF] = for e^> Sp , (2.9) 

where Sp is the Fermi energy. We assume that the temperature is zero. 

The Hartree-Fock potential (2.8) is seen to be state-dependent mainly due to the 
exchange term. Usually, the exchange term is replaced by a local density approximation 
[3, 24]. If we neglect the exchange term altogether and suppress the exclusion of the state 
'^K = y^u in the remaining Hartree potential, we obtain 

t/e(r)^.(f) = / Sr'P^P^ , (2.10) 



^We leave away an explicit notation of the spin degrees because spin-dependent interactions are 
neglected. 
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pf W) = E^«eo(/P^(rO ^^i?) . (2.11) 

K 

The solution of the remaining set of selfconsistent Hartree equations (Eqs. (2.5) with 
(2.6), (2.7), (2.10), (2.11)) is still a considerable technical problem. It has been carried 
through for instance in Ref. [3]. 

A considerable simplification is obtained if the selfconsistent average potential (2.6) 
is replaced by a phenomenological ansatz. Most of the shell structure calculations for 
compact metal clusters have been performed on the basis of simple phenomenological 
potentials, in particular the Nilsson potential [9], the Saxon- Woods potential [7, 11], 
and the ,, wine-bottle" potential [7]. It can indeed be seen from the results of Ref. [3] 
that for atom numbers n > 40 the selfonsistently calculated potential V{r) resembles 
a Saxon- Woods potential and roughly even to a square well. We, therefore, felt justified 
to represent the average potential V{r) for the case of a bubble cluster by an infinite 
square well with the boundaries given by the distribution of the positive ions 

nr) = |:^° ^^f^<'<R^ . (2.12) 
^ ' I -l-oo otherwise ^ ' 

This simple choice of the shell model potential has the great advantage that the eigen- 
functions of the Schrodinger equation (2.5) are linear combinations of spherical Bessel- 
and Neumann functions and the eigenenergies e^, are easily obtained from the boundary 
conditions at r = Ri^2 [23]. 

The well depth Vq in (2.12) is of the order of 0.5 Ry (1 Ry = 13,6 cV). Its value is 
not relevant for our results because the shell correction energy (see Eq. 2.16) turns out 
to be independent of the constant Vq. 

We now have to determine the total energy E of the metal cluster as a function of 
the variable / (see Eq. (2.4)). Given the fact that the ion density Pion(^) can be considered 
to be constant inside the matter distribution (see (2.1)), we may write the energy of the 
system as a sum of the energy E^y) of a ,, liquid drop" and a ,, shell correction energy" 
-E'sheii following Strutinsky [13] 

-E'tot — -E'ld + -E'sheii (2-13) 

The liquid drop energy can be written as a sum of a (negative) term proportional to the 
volume V of the system and a (positive) term proportional to the surface S. 

Since we consider (positively) charged clusters, we have to add the electrostatic 
energy Ech of the sytem 

Eld = -tY + aS + Ecy, (2.14) 

which turns out to be a much bigger term than in a neutral cluster. 

Some authors [10] include a term proportional to the average curvature of the surface 
which is the 3'''^ term in the expansion of the energy-density functional of a leptodermons 

sytem in terms of n~^^^. In the case of a bubble shape, the curvature terms arising from 
the inner and outer surface have opposite signs beside the fact that their absolute value is 
much smaller than the corresponding surface term. We omit the curvature term thereby 
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following Ref. [21]. Consideration of the curvature term would favour spherical bubbles 
as compared to compact spheres. 

The ,, macroscopic" electrostatic energy Ecb is given as a function of smooth density 
distributions of the positive and negative charge by 



Cb 



f[Pe{r) - Pion(r)][pe(r') - Pion(r')] 



(2.15) 



and the shell-correction energy E'sheii as a function of the single-particle energies e^, by 

-E'sheii = ~ ' (2-16) 

K 

where the occupation probabilities were defined in (2.9), whereas the quantities Uk 
represents smooth occupation probabilities in a cluster in which the shell structure is 
washed out using the Strutinsky prescription [13]. 

In our simple phenomenological model, we represent the distribution of the positive 
surplus charge [pion('^) — Pe('")] by a simple ansatz: From classical electrodynamics we 
know that the surplus charge ought to be localized at the outer surface of the metallic 
bubble. We, therefore, assume the total (positive) surplus charge qcq to be distributed in 
a thin layer of thickness e along the outer surface 



where 



[Pion(r) - pe(r)] = Sp ■ eo{r - (i?i - e)) ^o(i?i - e) , 

. 3geo 

" 47r[i?f - {R, - £)3] ■ 

Calculating the Coulomb energy (2.15) for this distribution we obtain 



(2.17) 



Ecb — 



2--4 

Ri Rl 



(2.18) 



Up to terms of the order of , this can be written 

Ech 



(f4 
2Ri 



1 e I 



2 



(2.19) 



The term of order (^-^ in (2.19) represents the Coulomb energy in the hmit that the sur- 
plus charge is located in an infinitely thin layer at the outer surface. One sees from (2.19) 
that the Coulomb energy increases if the surplus charge is distributed homogeneously in 
a layer of finite thickness e <^ Ri. Consequently, we used the Coulomb energy in zeroth 



order of in our calculations. 

The total energy of the spherical metallic bubble has thus the form 

47r 



-E'tot — 



-(i??-i?^) + a47r(i?? + i?^) + g 



+ Eshewif) ■ 



(2.20) 
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Subtracting from this expression the energy of a compact spherical cluster of the same 
charge qe^ and the same volume we obtain 

AE{f;q): = AEld(/; g) + AE,heii(/; n - g) , (2.21) 

where 

AEUf; Q) = 47ra (Rj + Rl - Rl) + (^1- - i-) (2.21') 

and 

A£;,heii(/; n - 1) = £;sheii(/; n-q)- £;sheii(0; n-q) . (2.21") 

We note that the single particle potential (2.12) becomes a simple central square- well 
in the limit R2 ^ (i.e. / ^ 0). Thus the shell correction energy -Eshcii(/ = 0) is 
obtained by substituting the eigenvalues £«(/ = 0) in a simple central well with infinite 
wall at r = i?o- The occupation probabilities have to be chosen in each case according to 
Strutinsky's prescription. 

The 1®* term and 2^^ term of (2.21') are simple functions of the parameter / due to 
the relations 

/ 1 xl/3 

Ri-^o- [jZrj) ' (2.22) 

R2^Ro- (737) ^ • (2-22') 
The radius Rq of the compact spherical cluster is related to the number of n of atoms by 

n^/^ (2.23) 

with the radius parameter given in (2.2). 

The surface tension a and the radius of the Wigner-Seitz cell are thus the only 
parameters which specify a given metal in our model. 

The difference A£'ld(/; q) between the LD-energies of the compact cluster and the 
bubble (Eqn. (2.21')) is seen to consist of a positive term describing the increase of 
the surface energy and a negative term which represents the reduction of the repulsive 
Coulomb energy. As far as the ,, macroscopic" part of the energy is concerned, a preference 
for the bubble geometry may only occur through the reduction of the Coulomb energy 
which increases with the surplus charge q of the metal cluster. For a given metal and 
a given ,,size" n of the cluster there will thus be a critical value qc of the charge at which 
the (spherical) bubble becomes the configuration of lower energy. It will depend on the 
type of the metal (i.e. on the value of the surface tension a and the radius parameterr^) 
whether the cluster is still stable versus fission or emission of atoms at this value of the 
charge. 

The difference Ai^sheii of the shell correction terms can be positive or negative. If 
the number (n — q) of valence electrons happens to be a magic number for the bubble 
geometry, and not for the compound spherical form, AE'gheii will be a negative number 
and thus favour the formation of a bubble a.v.v. 
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The question whether, for given charge, the bubble has a lower energy than the 
compact sphere, depends sensitively on the value of the surface constant a and the radius 
parameter r^. More precisely, it depends on the "fissility" parameter which is defined to 
be the ratio of the Coulomb energy Eq\, and twice the surface energy Es- The factor 2 is 
inserted in order to retain the definition used in nuclear physics: 

V _ Ec^U) _ q'el _ ( 1-f \ 

^ - 2Es{f) - 16na ■ R,{Rl + " Al + f^' ) ' ^ ^ ^ 

with Xq, the fissility parameter of a compact sphere, being defined by 

Ec^if = 0) q'el 
^0 = TTETTT ^ = • (2.25) 



2Es{f = 0) i67ra r 



n 



The difference Aii^LD between the LD-energy of the spherical bubble and of the compact 
spherical cluster (see (2.21')) measured in units of 2Es{f — 0) is given by the function 



^ ^AE^j,if;q) _ 1 



2Es{f = 0) 2(1 - /)V3 



1 + f /3 - (1 - /)2/3l _ Xo [l - (1 - /)i/3l _ ^2.26) 



of the hole parameter / (0 < / < 1), the term on the r.h.s. is positive and the 2^'^ 
term is negative. Clearly, for large enough Xq, the function F becomes negative, i.e. the 
bubble shape corresponds to a lower energy. The stationarity condition 

dF 

Qj - (2.27) 

has the explicit form 

l + f/'-XQf/^{l-f) = 0. (2.28) 

Solving Eq. (2.28) one can find that the metallic cluster with charge q ■ Cq and the atom 
number n have their lowest energy for a spherical bubble shape if the fissility parameter 
Xo < 3,4 . 

A crucial question is the dependence of the total energy on deformation: It can 
be easily shown that the LD-energy of a bubble decreases as a function of deformation 
because the decrease of the repulsive electrostatic energy turns out to be greater than 
the increase of the surface energy. A spherical bubble may, however, be stabilized against 
deformation by shell effects. If the {n — q) valence electrons of the charged metal bubble 
correspond to a closed shell configuration in the spherically symmetric potential (2.12), 
the shell energy yields an additional binding of a couple of eV. As one deforms the bubble 
the absolute value of this negative shell energy decreases as a function of deformation. 
In this way a barrier against fission is produced. The calculation of the total energy of 
the bubble as a function of the deformation implies that we determine the eigenvalues 
in a deformed bubble potential. This is a difficult task if we were to tackle it in full 
generality. There is, however, a family of deformed shapes which can be transformed into 
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a spherical shape by a scahng transformation. For this "scahng model", it is relatively 
simple to calculate the eigenvalues: 

Assume that the outer (5*1) and inner (5*2) surface of the deformed bubble are 
concentric spheroids with the half-axes 01(2) and ci(2) 

Sm-. ^^ + -^ = 1 (2.29) 

"1(2) 4(2) 

enclosing the constant volume ^(-Ri — R2) 

Rl^Rl-Rl^ a\ ci - a\c2 . (2.30) 
An infinite square well with the boundaries (2.29) 

t>(x,y,^) = -yo 770(f) , (2.31) 

where 



?7o(f) 



1 for X G volume Q enclosed by S\ and 5*2 
—00 otherwise 



transforms into a simple spherically symmetry potential of the type (2.12) which depends 
on the scaled variable 

i = \x- ri = Xy; ( = l^z (2.32) 

or the corresponding polar variables p, 7?, </? 

^ = p sin 1? cos </?; 77 = p sin sin y?; ^ — pcos'&. (2.32') 

The 2 scaling parameters A, p are related to each other by the constraint (2.25) which 
takes the form 

AV = 1 • (2.33) 

Thus they can be expressed by a single deformation parameter. Using the deformation 
parameter 5 introduced by S. G. Nilsson in Ref. [25] the scaling parameters are given as 
a function of 5 by 

In the scaled variables, the Schrodinger equation takes the form 

{Ho + Hi)^u = e,ip, . (2.35) 
The hamiltonian Hq is spherically symmetric in the ^, 77, C coordinates 

Ho^-^^\r,C + V{p) , (2.36) 
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where A^^^ is the Laplacian operator in ^, r^, C space and 



Ms ^M{1 + '^6)1(1 



(2.37) 



The scaled square well potential is now spherical 

V{P) = 



-Vo for R2 < p < Ri 
+00 otherwise 



(2.38) 



The term Hi in (2.31) represents the deformation dependent part of the hamiltonian 

92 \ 



/7 - - — 



(2.39) 



For determining of the stability of spherical bubble with respect to elongation, we 
only need to consider small deformations 5. Consequently, we may treat Hi as a pertur- 
bation. This means that the single-particle energy can be approximately obtained as 
a function of the eigenenergies and eigenfunctions (^^ of the Hamiltonian Hq as 



+{^^ \Hi\ 



(2.40) 



For larger deformations, one has to diagonalize Hi in the basis of the s.p. states (p^,. In 
all the results shown in Section 3, the eigenvalues e^, were determined by diagonalization 
of Hi in a sufficiently large subspace of s.p. states </?^. 

The matrix-elements of Hi can be easy evaluated using the following artifice 



92 



92 



1 



1 f^Ms 



dC di"^ dv"^ 8 
where is the Lagrange polynomial. Simple algebra leads to the expression 



{^u \Hi\ 



6 /2M, 



{eu - e^f{^^ \p^P2{cose) \ , 



(2.41) 



(2.42) 



The matrix elements of where evaluated numerically, while the matrix elements of 
spherical harmonics are expressed in terms of Clebsch-Gordan coefficients. 



3 Results 

The calculation was performed for metallic agglomerates of sodium (Na) and of antimo- 
nium (Sb). The liquid drop parameters are the same as in Ref. [21] where the fission of 
charged clusters was discussed. For Na clusters we have used the set 

2.070 A , T = 0.03017 eV , a = 0.01894 eV 
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and for Sb clusters 

r,= 1.130 A , T = 0.4552 eV , a = 0.02474 eV . 

The results for the charged sodium clusters are presented in Fig. 1 The number n 
of atoms in clusters vary from 700 to 10000, while the charge number q changes from 30 
to 200. The liquid drop estimate of the binding energy {Eld, Eq- 2.1), the energy gain 
{AElo, Eq. 2.21') with respect to the energy of a compact spherical cluster, the fissility 
parameter (X, Eq. 2.24) and the equilibrium hole fraction (/, Eq. 2.4) are drawn. It is 
seen in the Fig. 1 that only massive and relatively highly charged Na clusters favour the 
bubble solution. The average energy gain is around 0.3 eV per atom. The hole fraction 
parameter / varies from 0.3 to almost 0.9 and the hole in the bubble is larger the higher 
is the charge of the cluster. Unfortunatelly the fissility parameter X is always larger than 
one. This means that there is no fission barrier in the liqud model and there is a little 
hope that the shell effect is large enough to stabilize the sodium bubble. 

A more optimistic situation exists for the Sb clusters, which have a smaller Wigner- 
Seitz radius constant. In Fig. 2 we plotted analogous results as in Fig. 1. One can see 
that already small and not highly charged clusters prefer the bubble configuration. The 
energy gain with respect to the compact shape could even reach 1 eV per atom, also 
the binding energies per atom are few times larger than for Na bubbles. But the fissility 
parameter X is here also larger than one. One can conclude from the LD results presented 
in Figs 1 and 2 that within the LD model stable bubble clusters are rather improbable. 
Only the shell effects may stabilize them with regard to fission. 

The electronic scheme obtained for the infinite spherical square well potential (2.12) 
is plotted in Fig. 3 as a function of the hole fraction / (2.4). The energy unit used takes 
into account that the eigenvalues of the infinite square well scale with •nr'^l'^ . The radius 
constant for Sb is used here. It is seen in the figure that for larger / the orbitals with 
the node number n > 1 corresponds to much higher energies than those with n = 1 and 
the levels are well separated from the others. This effect leads to strong shell effects. A 
similar tendency was also observed in the shifted harmonic oscillator [23]. New magic 
numbers corresponding to the orbitals with l—O, 1, 2, 3, .... are found: 

Mbubb. = 2, 8, 18, 32, 50, 72, 98, 128, 162, 200, 242, 288, 338, 392, 450, 512, 578, 648.... 

These magic numbers are quite different from those observed in the compact spherical 
clusters [26] : 

Mcomp. = 2, 8, 20, 34, 58, 92, 138, 168, 254, 338, 438, 440, 542, 556, 676 

The Strutinsky shell correction energy for the Sb bubble cluster with the hole fraction 
f—0.7 is plotted as a function of n in Fig. 4. The shell correction energy for a magic 
number of valence electrons in the bubble cluster is negative and reaches even -3 eV for 
n=648. The magnitude of the shell correction is comparable with that for supershells for 
heavy clusters with 1074 < n < 3028 discussed by Brack in Ref. [26]. Such a huge shell 
effect for magic clusters should protect them against fission. It can be seen in Fig. 5, 
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where the fission barrier for the magic bubble cluster 2ioSb'^^'^ is analysed. In the upper 
l.h.s. part of Fig. 5, the liquid drop energy is drawn as a function of the hole fraction /. A 
pronounced bubble minimum is observed at / ^ 0.4 . The liquid drop part of the energy 
at /=0.4 is plotted as a function of deformation parameter 6 in the upper r.h.s. part of 
Fig. 5. As one can expect there is no fission barrier in this case. The two lower graphs 
represent the shell energy (2.16) and the fission barrier beeing the sum of the LD part and 
the shell correction (2.13) as a function of S. It is seen that Esheii=-'^-i'^ eV for spherical 
configuration and its magnitude decreases with deformation and oscillates. This effect 
produces a prononced fission barrier which makes fission of 2iO'S'&'''^° rather improbable. 

4 Summary and discussion 

We have investigated the shell correction energy of charged spherical bubble clusters. The 
lower limit of mass and charge numbers, where bubble clusters begin to exist depends 
sensitively on the choice of the LD parameters. The investigation shows that promising 
candidates for bubble structure are Sh clusters with a magic number of itinerant electrons 
in the charged agglomerate. The charge number is found to be a fraction of 0.04 to 0.15 
of the atom number. The Na cluster should contain large atom number, preferably a few 
thousands and they should be highly charged {q fa 100). 

We have found strong shell effects which may give rise to shell energies of up to -4 eV 
for certain magic numbers. By calculating the LD-energy for deformed bubbles and the 
deformation dependence of the shell effect, we found that the fission barriers are of the 
same order of magnitude as the shell-correction energy for the spherical bubble solution. 
In favorable cases of magic numbers, this is sufficient to reduce significantly the probability 
for spontaneous fission. The origin of the shell effects is the high degeneracy of orbitals 
with large angular momentum and the rapid energetic increase of radial oscillation modes 
as a function of increasing bubble radius. 

We expect that similarly big shell effects could exist in the clusters in which the 
outer layer from other material is added. If would exist the potential barrier between 
the inner and outer material then the electrons in the outer metallic layer could feel the 
magic numbers like those found for the bubble cluster. It would be a great challenge for 
experimentalists to produce (by epitaxy?) such objects. 
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Figures captions: 



1. Liquid drop estimate of the binding energy (-Eldj Eq. 2.1), energy gain {AEld, Eq. 
2.21') with respect to the energy of a compact spherical cluster, fissility parameter 
(X, Eq. 2.24) and the equihbrium hole fraction (/, Eq. 2.4) for the charged {q ■ Cq) 
sodium cluster with the bubble structure as function of n. 

2. The same as in Fig. 1 but for charged clusters of antimonium. 

3. Electronic level scheme as a function of the hole fraction / = {R2/R1Y for the 
infinite spherical square well (2.12). The energy unit used takes into account that 
the eigenvalues of the infinite square well scale with The radius constant for 
Sh cluster is used here. 

4. Shell correction energy for the Sh bubble cluster with the hole fraction /=0.7 as a 
function of n. 

5. Liquid drop energy of 2wSh'^^^ as a function of the hole fraction / (upper l.h.s. fig- 
ure) and the deformation parameter 5 (upper r.h.s). The two lower graphs represent 
the shell energy (2.16) and the fission barrier (2.13) as a function of 5. 
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Figure 3: 
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